A numerical solution to the local cohomology problem in U(l) chiral gauge theories 
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We consider a numerical method to solve the local cohomology problem related to the gauge 
anomaly cancellation in U(l) chiral gauge theories. In the cohomological analysis of the chiral 
anomaly, it is required to carry out the differentiation and the integration of the anomaly with 
respect to the continuous parameter for the interpolation of the admissible gauge fields. In our 
numerical approach, the differentiation is evaluated explicitly through the rational approximation 
of the overlap Dirac operator with Zolotarev optimization. The integration is performed with 
a Gaussian Quadrature formula, which turns out to show rather good convergence. The Poincare 
lemma is reformulated for the finite lattice and is implemented numerically. We compute the current 
associated with the cohomologically trivial part of the chiral anomaly in two-dimensions and check 
its locality properties. 

PACS numbers: 11.15.Ha, 12.38.Gc. 



I. INTRODUCTION 

The construction of gauge-covariant and local lat- 
tice Dirac operators which satisfy the Ginsparg- Wilson 
relation[aiJ|iaail3, 



j 5 D + D l5 = 2aD 75 D, 



(1) 



has made it possible to introduce Weyl fermions on the 
lattice and construct anomaly-free chiral gauge theories 
with exact gauge invariancepl 1^. TToL ITU [l3| . One of 
the crucial steps in the gauge-invariant construction is to 
establish the exact cancellation of the gauge anomaly at 
a finite lattice spacing. 

In the case of U(l) chiral gauge theories^ |37j. the ex- 
act cancellation has been achieved through the cohomo- 
logical classification of the chiral anomaly [cl l35l l3f| [3I I . 
The anomaly is given in terms of lattice Dirac operator |3 

El mill El Ea as 



q(x) = tr{7 5 (l - aD){x,x)} 



(2) 



and the local field q(x) is a topological field in the sense 
that it satisfies 



8q{x) = 



(3) 



under a local variation of the gauge field. It follows from 
this property that the anomaly is cohomologically trivial, 



Here k^(x) is a gauge-invariant local current. This local 
current is in turn used in the gauge-invariant construction 
of the functional measure of the Weyl fermions. 

If one thinks of the practical computation of observ- 
ablcs in the lattice U(l) chiral gauge theories, it is re- 
quired to compute the local current in Eq. for every 
admissible gauge field. The purpose of this paper is to at- 
tempt a numerical computation of the local current. We 
will compute k^(x) in two-dimensions numerically and 
check its locality properties. We will also check how the 
exact cancellation of gauge anomaly works on the finite 
lattice. 

This paper is organized as follows. In section [H] we 
formulate the vector-potential-representation of the link 
variables for the admissible U(l) gauge fields on the fi- 
nite lattice. Using this representation, we introduce one- 
parameter families of the admissible fields for the in- 
terpolations. In section IIIII we describe our numerical 
method to compute the bi-local current which is the first- 
differential of the chiral anomaly with respect to the vec- 
tor potential. In section Hvl the Poincare lemma is refor- 
mulated for a finite lattice so that we can carry out the 
cohomological analysis directly on the finite lattice. The 
result of the cohomological analysis in two-dimensions 
is summarized. In section [V] we describe our numerical 
result of the computation of the local current (x) . Sec- 
tion is devoted to a summary and discussions. 



52e a f(x) = &lk M (x), q a {x)=q{x)\ u ^ uea , (4) 

a 

for an anomaly- free multiplet of Weyl fermions which sat- 
isfies the condition of the U(l) charges, 

a 
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II. ADMISSIBLE U(l) GAUGE FIELDS 
ON A FINITE LATTICE 

Our first step is to formulate the vector-potential- 
representation of the link variables associated with an 
admissible U(l) gauge field on a finite lattice. Such a 
representation has been formulated in the original co- 
homological analysis in However, it is constructed 
for the admissible gauge fields on the infinite lattice and 
the resulted vector-potentials are not bounded in general. 
This representation therefore does not seem to be useful 



2 



for numerical implementations. But, as has been shown 
in our previous paper jS^, it is possible to formulate the 
bounded and periodic vector-potential representation for 
the admissible gauge fields on the finite lattice. 

We set the lattice spacing a to unity and consider U(l) 
gauge fields on a finite two- or four-dimensional lattice 
(n = 2,4) of size L with periodic boundary conditions. 
L is assumed to be an even integer for simplicity. The 
gauge fields on such a lattice can be represented through 
periodic link fields on the infinite lattice, 

%/i)eU(i), x = {x 1 ,x 2 ,--- ,x n ) ez n , (6) 

U(x + Li), /i) = U(x, fi) for all fx, v = 1, • • • , n.(7) 

The independent degrees of freedom are then the link 
variables at the points in the region 

T n = {xeZ n \-L/2<x^<L/2}. (8) 

As to gauge transformations 

U(x,n) -> A(x)U(x,li)A(x + fi)- 1 , (9) 

we consider only periodic functions A(ir) S U(l) which 
preserves the periodicity of the link field. 

We impose the admissibility condition on the U(l) 
gauge fields: 

\Ffj,v(x)\ < e for all x,(jl,v, (10) 
where the field tensor F /lu (x) is defined through 

F^( X ) = ihlP^(.T), -7T < F MV (X) < 7T, (11) 

p^A x ) = u(x,(i)u(x + ji,i')u(x + v,iiy 1 u(x,vy 1 . 

(12) 

We require this condition because it ensures that the 
overlap Dirac operator 0,3, which we adopt in this work, 
is a smooth and local function of the gauge field for 
|1— P liV (x)\ < 1/30 in four-dimensions and 1 1 — P fll ,(x) \ < 
1/5 in two-dimensions p| 56]. For e < 7r/3 the admissible 
U(l) gauge fields on the finite lattice can be classified 
uniquely by the magnetic fluxes to m „ (integers indepen- 
dent of x) where 

1 L_1 

"V, = ^— F fll ,(x + sp, + tv). (13) 

s,t=0 

In this respect, the following field is periodic and can be 
shown to have constant field tensor equal to 27rm M „/ 1 L 2 : 

V [m] (x, fi) = e -f¥[ il5 ^^-i E,>m ™^+Y,„<„ i 

(14) 

where the abbreviation i M = x^ mod L has been used. 
Then any admissible U(l) gauge field in the topological 
sector with the magnetic flux m^ v may be expressed as 

U(x,/J,) = U(x,n) V[ m ](x,fj,). (15) 



We may regard U (x, j£) as the actual local and dynamical 
degrees of freedom in the given topological sector. This is 
because the magnetic flux m^ u is invariant with respect 
to a local variation of the link field. 

The following lemma shows that it is possible to es- 
tablish the one-to-one correspondence between U(x, /i) 
and a periodic vector potential with the desired locality 
properties on the finite lattice |55|. 

Lemma [H] There exists a periodic vector potential 
A fl (x) such that 

e lA ^ x) = U{x,n), (16) 

dMx) - d„A„(x) = F^(x) - — jf-, (17) 

f L4 M (»|< 7r(H-4||a: ||) \\x\\<Lj2 
\ I^WI < n ( l + 2L + 2 ( n - l ) L2 ) otherwise. 

(18) 

Moreover, if A'^{x) is any other field with these proper- 
ties, we have 

A' ll (x)=A ll {x)+d^{x), (19) 

where the gauge function lj(x) takes values that are in- 
teger multiples of 2ir. The explicit formula of A^{x) in 
two-dimensions is given in the appendix. 

As emphasized in 0, it is important to note that the 
locality properties of gauge invariant fields should be the 
same independently of whether they are considered to 
be functions of the link variables or the vector potential. 
Since the mapping 

A^x) -> U(x,fi) = e iA »W (20) 

is manifestly local, this is immediately clear if one starts 
with a field composed from the link variables. In the 
other direction, one may start from a gauge invariant 
local field 4>(y) depending the vector potential. Then the 
key observation is that one is free to impose a complete 
axial gauge taking the point y as the origin. Around y 
the vector potential is locally constructed from the given 
link field and 4>(y) thus maps to a local function of the 
link variables residing there. 

The vector potential A fJ- (x) represents an admissible 
field through e lA ^ x "> x V[ m ](x,/J,) and the associated field 

tensor F tiv (x) = d tl A v (x) — d^A^x) + 2 "" A "' is hence 
bounded by e. It is straightforward to check that this 
property is preserved if the potential is scaled by a factor 
t in the range < t < 1, i.e. we can contract the vector 
potential to zero without leaving the space of admissible 
fields. Then, in the cohomological analysis of the chiral 
anomaly, we may choose Vj m i (x, fi) as the reference gauge 
field in a given magnetic flux sector and may consider 
the interpolation between the arbitrary gauge field in the 
same sector and the reference field as follows: 

U t (x,fi) = e ltA ^ x) x V [m] {x,ii) t G [0,1]. (21) 
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III. CHIRAL ANOMALY AND ITS 
TOPOLOGICAL PROPERTIES 



the bi- local current, we introduce the parameter repre- 
sentation of the inverse square root of as 



For our numerical application, we adopt the overlap 
Dirac operator[2, H given by 



d = - 2 1 



H„ 



75" 



where H w is the Hermitian Wilson-Dirac operator, 



mo 



(22) 



(23) 



(0 < too < 2). It is a smooth and local function 
of the admissible gauge field for |1 — P liu {x)\ < 1/30 
in four-dimensions and for |1 — P liv {x)\ < 1/5 in two- 
dimensions 0(5^. Then the chiral anomaly, Eq. J2J), is 
given in terms of H w as 



q{x) = --tr 



(x,x) 



(24) 



As is clear from this expression, q(x) is a topological field 
with respect to any local variation of the admissible gauge 
field, 



0. 



(25) 



This topological property of the chiral anomaly can 
be cast into properties of a certain gauge-invariant bi- 
local current. This bi-local current is defined through 
the differentiation and integration with respect to the 
continuous parameter t for the interpolation as follows: 



n \dA v {y) 



(26) 



A~>tA 



The original topological field can be expressed with the 
bi-local current as 

q(x) = q[m]{x) + ^2 3A x ,y)Av{y), (27) 

where q[ m ](x) denotes the topological field for V[ m ](x,ij). 
Then the topological property and the gauge-invariance 
of q{x) imply that j v {x,y) satisfies 



(28) 



xer„ 



These conditions provide the initial conditions for the 
cohomological analysis. 

For our purpose, it is required to compute the above bi- 
local current numerically keeping the conditions Eq. Ij28(l 
within a given accuracy. Our strategy is the following. 
First of all, in order to obtain more explicit formula of 



dt H w 
TTt 2 + HI 



(29) 



Then we can perform the differentiation of q(x) with re- 
spect to the vector potential A^(x) explicitly |5q| as 



dq{x) 
dA v (y) 



1 

"2 tr 



dt 



1 



7T i 2 + HI 



{t 2 V^y) - H w V^(y)H M 



1 



t 2 + HI 



[x,x) 



(30) 



where 



1 



V^(y) = -im-Hl-i^Sx^x+^zU^x) 

+ (1 + j^Sy^Sx^+^U^z)- 1 } . 



(31) 



For the numerical evaluation of the differentiation, 
we then adopt the rational approximation of the 
overlap Dirac operatorJH^ E3] with the Zolotarev 
optimization|6ll l62j|. Namely, q(x) is approximated by 



the formula with the degree N r 



1 f Nr 

{ k=l 



bk 



hi 



C2fc-1 



-(x,x) 



(32) 



h w is defined by H w /X mm where A m i n is the square root 
of the minimum of the eigenvalues of H 2 . The definitions 
of the coefficients Cfc and bk are given in the appendix. 
With this approximation, Eq. (|30|1 can be approximated 
as follows: 



dq(x) 
dA v {y) 



1 



-tr 



1 



h 2 



C2fc-1 



{c2k-iVfj,{y) - h w v f _ l (y)h w ) 



1 



hi + C 2 k-1 



(x,x) 



(33) 



where v^(y) = V^(y)/A m i n . Finally, for the integration 
back with respect to the parameter for the interpolation 
in Eq. l|26|l . we adopt the Gaussian Quadrature (Gauss- 
Legendre) formula with N g points: 



dq(x) 
dA u {y) 



(34) 



A^UA 



where {(U, Wi)\i = 1, • • • , N g } is the set of the abscissas 
and weights of the degree N g . 

The formulae, Eqs. I|33|l and (|34|) . are still complicated 
for the numerical computation. But, in two-dimensions, 
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it is not demanding numerically for relatively small lat- 
tice sizes to diagonalize H w and store all the eigenvectors. 
Namely, we can write 

h w (x, y) = J2 M x ) (tM (35) 

and then it is possible to evaluate Eqs. explic- 
itly. 

In this approximation, the gauge-invariance of the bi- 
local current and the second property of Eq. (|2*5)l can be 
preserved exactly. As to the first property of Eq. (|28|l . 
we will see below that the choice N r — 18 and N g = 20 
gives good convergences for the admissible gauge fields 
on the two-dimensional lattice of the size L — 8, 10, 12 
and the original topological field, q(x), can be reproduced 
through Eq. I|27|) within the error less than 1 x 10~ 14 . 

IV. COHOMOLOGICAL ANALYSIS OF 
CHIRAL ANOMALY ON A FINITE LATTICE 

A. The modified Poincare lemma on a finite lattice 

Our numerical cohomological analysis of the chiral 
anomaly will be performed directly on a finite lattice. 
For this purpose it is required to reformulate the Poincare 
lemma on the lattice 8] for a finite lattice. For the de- 
tail of the proof of the lemmas, refer to our previous 
paper[55j. 

In the following we will consider tensor fields f^ - ^ (x) 
on r„ that are totally anti-symmetric in the indices 
111,-" , /ife- Such tensor fields may be regarded as pe- 
riodic tensor fields on the infinite lattice, 

fm-n k ( x + L fy = Ui-Hk( x ) f OT all fi,v = 1,- ■ ■ ,n. 

(36) 

The locality properties of such fields are assumed to be 
as follows: for a certain reference point xq 6 T n and 
-L/2 < (x„ - x 0fl ) < L/2 (mod L), 

< da+lli-ioDe-ll*-^, 

(|| a; -io ||<£/2), (37) 

(\\x-x \\>L/2). (38) 

g is a localization range of the tensor field, d and pi > 
are certain constants that do not depend on L. || x — xq || 
is the taxi driver distance from xo to x. This local- 
ity properties hold true for the differentials of the chiral 
anomaly given in terms of overlap Dirac operator 3 
with respect to the admissible gauge fields [q [o3|. 

The differential forms on the finite lattice are intro- 
duced as in the continuum, following || . If we adopt the 
Einstein summation convention for tensor indices, the 
general fc-form on r„ is then given by 

/ 0) = tt/mi-wc {x)dx^ ■ ■ ■ dx^ k . (39) 



The linear space of all these forms is denoted by fife. An 
exterior difference operator d : — > flk+i may now be 
defined through 

df{x) = — { d IJ ,f lll ,...,» k {x)dXi l dx lll ■ ■ ■ dx^ , (40) 

where 9 M denotes the forward nearest-neighbor difference 
operator. The associated divergence operator d* : ilk — » 
flk-i is defined in the obvious way by setting d* f = if 
/ is a 0-form and 

d*f(x) = 1 ^ d*f^ 3 ... lih (x)dx li3 - --dx^ (41) 

in all other cases, where d* is the backward nearest- 
neighbor difference operator. 

By definition, the divergence operator satisfies that 
d* 2 = and therefore the difference equation d* f = is 
solved by all forms / = d*g. It has been shown that in the 
infinite lattice these are in fact all solutions, an exception 
being the 0- forms where one has a one-dimensional space 
of further solutions^. This result is the lattice counter 
part of the Poincare lemma known in the continuum the- 
ory. On the finite periodic lattice the lemma does not 
hold true any more, because the lattice is a n-dimensional 
torus and its cohomology group is now non-trivial. How- 
ever, the lemma can be reformulated so that it holds 
true up to exponentially small correction terms of order 
0(e~ L/2e ) and for the form satisfying J2x^f( x ) = °> 
it holds exactly even on the finite lattice. [63 The lat- 
ter result is the lattice counter part of the corollary of 
dc Rham theorem known in the continuum theory. The 
precise statements are the following. 65J 

Lemma HVl a (Modified Poincare lemma) 

Let / be a fc-form which satisfies 

d*f = and f( x ) = if k = °- ( 42 ) 

Then there exist a form g € £lk+i and a form A/ G Qk 
such that 

f = d*g + Af, \Af^... n (x)\ <cL°e~ L ^. (43) 

Lemma HVl b (Corollary of de Rham theorem) 

Let / be a fc-form which satisfies 

d*f = and f( x ) = °- ( 44 ) 

Then there exist a form g £ Qk+i such that 

/ = d*g. (45) 

The explicit formula of 3(2;) in two-dimensions (n = 2) 
is given in the appendix. 
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B. A solution to the local cohomology problem 
on the two-dimensional finite lattice 

In two dimensions, the cohomological analysis results 
in the following formula for q(x): 

q(x) = A(x) + d*h„(x) + Aq(x), (46) 

where 

A{x) = q [m] {x) + (f>^(x)F^(x), d*^ v (x) = 0. (47) 

<Ptj,v(x) is obtained through the applications of the mod- 
ified Poincare lemma and the corollary of the de Rham 
theorem: 

j v {x,y) ^ e vtl (x,y)d*+Aj v (x,y), (48) 

^A>(z,x) ¥± -2A^(x)^ (49) 
zer 2 

and 

<j>n V (z) = i ^2 0^(z,x) - AQ^x). (50) 

zer 2 

In two-dimensions, 4>^, v {x) is a constant and assumes the 
form ^^{x) — 7[ m , u) ]e f ji/. h fJ _(x) is obtained through the 
application of the modified Poincare lemma: 

9^(x,y) - S x . y ^2 9^(z,y) d\T Xllv {x,y) 
zGr 2 

(51) 

and 

Ki x ) = \ ^Z T ^ P { x ^y)Fyp{y)- ( 52 ) 

yer 2 

Aq(x) is defined by 

Aq(x) = A$ llv {x)F^(x)+^Aj v {x,y)A v {y). 

(53) 

For the anomaly-free multiplet satisfying the condi- 
tion J2 a e a — 0, the anomaly part A(x) — m m i (x) + 
l[m\^^F lll j(x) cancels up to an exponentially small term, 

J2 e *A a (x) = AA(x), (54) 

a 

where A a (x) = A{x)\ u ^ Ueot . Since £x £ r 2 AA(x) 
. we may apply the modified Poincare lemma IV. a to 
obtain 

AA{x) ^ d^Ak^x). (55) 

On the other hand, since XLer 2 ^<?( x ) = 0' we may also 
apply the lemma IV. b as follows: 

Aq(x) — d;Ah^x). (56) 



Then we can see that q(x) is cohomologically trivial, 

J2e a q a (x)=d;k,(x), (57) 

a 

where k^(x) is given explicitly as 

k M {x) = J2 e <x{K{x) + Ah„(x)} a + Ak^x). (58) 

a 

Therefore, once the bi-local current jn(x,y) is com- 
puted, the current fc M (x) can be obtained by the sequence 
of the applications of the lemmas IV. a, b. The numeri- 
cal implementation of this step is straightforward using 
the explicit solutions, g(x), of the lemmas given in the 
appendix. 

V. NUMERICAL RESULTS 

We now describe our result of numerical computations 
of the local current fc M (x) . We consider the lattice sizes 
L = 8, 10, 12. Admissible gauge fields are generated by 
Monte Carlo simulation using the action 

g- 1 V F ^( x f -flV F a (kg) 

As reported in the topological charge is preserved 
during the Monte Carlo updates with this type of ac- 
tion, even when e is set to it. We adopt this option and 
check the locality of the topological field numerically for 
several values of f3. We consider the topological sectors 
with m 12 = 0,1 and the initial configuration is chosen as 
V[ m ](x,fx) with a given myi- 

For each admissible gauge field, we compute the vector 
potential A tl (x) formulated in section [HI We also com- 
pute the abscissas and weights {(ti,Wi)\ i — 1, • • • , N g } 
for Gaussian Quadrature formula with the degree N g . 
By this, the discrete interpolation of the given admissi- 
ble gauge field is fixed: 

{uW(x, fJ ,)=e it * A »W xV [m] (x,n) | i = l,-~,N g }. 

(60) 

In table m the abscissas and weights for the degree N g = 
20 are shown. 

For each (x, /i), all eigenvalues and eigenvectors 
of H w are computed numerically using the Householder 
method. We choose mo = 0.9 in H w . Then we compute 
the bi-local current j u {x 1 y) through Eqs. and lEHI) . 

We can check the convergence of the above procedure 
through Eq. I|27|) by computing the deviation 

X = max xe r 2 \S(x)\, (61) 

where 

5(x) = q(x) - q[ m ](x) - j v {x,y)A v (y). (62) 
yer 2 
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TABLE I: The abscissas and weights for the Gauss-Legendre 
formula are shown for N g — 20. t £ [0, 1]. 





U 


Wi 


1 


3.435700407452558E-003 


8.807003569576136E- 


003 


2 




2.030071490019346E- 


002 


3 


A 000070^07(100 TAO T7 1 AHO 

4.3882785874337U3.h-(J(J2 


3.133602416705449E- 


002 


4 


o.U441ol4UoooyUDlil/-UU2 


4.163837078835234E- 


002 


5 


0.120834046769925 


5.096505990862024E- 


002 


6 


0.181973159636742 


5.909726598075923E- 


002 


/ 


U.2440DD4yyU24ooD 


6.584431922458825E- 


002 


8 


c\ 0101 a r*c\ tt c r* j oo c\r\ 
U.olol4oy55o422y(J 


7.104805465919108E- 


002 


n 
9 


U.ooollHlH442yi ( ( 


7.458649323630189E- 


002 


1U 


U.4ol (oo( 3y43o2ol 


7.637669356536294E- 


002 


11 


0.538263260566749 


7.637669356536294E- 


002 


12 


0.613892925570823 


7.458649323630189E- 


002 


1 Q 

13 


0.686853044,50 II 10 


7.104805465919108-h- 


002 


14 


0.755433500975414 


6.584431922458825E- 


002 


15 


0.818026840363258 


5.909726598075923E- 


002 


16 


0.873165953230075 


5.096505990862024E- 


002 


17 


0.919558485911109 


4.163837078835234E- 


002 


18 


0.956117214125663 


3.133602416705449E- 


002 


19 


0.981985963638957 


2.030071490019346E- 


002 


20 


0.996564299592547 


8.807003569576136E- 


003 



Here the original topological fields q(x) and q[ m ](x) are 
constructed by computing all eigenvalues and eigenvec- 
tors of H w for the given admissible gauge held C/(x,/x) 
and for the reference gauge held V[ m ](x,fi), respectively. 
In this computation, we found that the topological holds 
have typically the values of order O(10~ 2 ) — O(10 -3 ). 
The integer topological charge Q = X^er 2 = 777-12 
is reproduced within the error of order O(10 -12 ) — 
O(10 -13 ). In table ITT1 we show the dependence of \ on 
the degree N g with N r = 18, L = 8 and (3 = 3.0 hxed. 

TABLE II: The convergence of the bi- local current is shown 
for various values of N g with N r = 18, L = 8, ft = 3.0 fixed. 

X 

4 O(10" 8 ) 
6 0(hn 12 ) 
8 O(10" 14 ) 

io o(hn 15 ) 

12 O(10 _15 ) 
16 O(10 _15 ) 



For the larger lattice sizes L = 10, 12, we found that % 
is less than 1 x 10~ 14 with the choice of N g = 20 and 
N r = 18. 

Once the bi-local current j^{x,y) is computed, the co- 
homological analysis of the chiral anomaly can be per- 
formed numerically by the sequence of the applications of 



the lemmas. The result can be expressed as (cf. Eq. (|46|) 

) 

q{x) = A{x) + d*h„(x) + Aq(x), (63) 

where 

A(x) (x). (64) 

In order to check the locality properties of the helds, 
q(x), A(x), h^{x) and Aq(x), we apply a small local vari- 
ation to the given admissible gauge held as U(x, fi) — > 
e l ^A x >U{x, n) where 

rj^x) = 0.05 x 27t4, Xo 5 Mi i (65) 

and compute the variations of the helds. For each varia- 
tion of the helds, S v f(x), we dehne 

8r)f(r) = max {\S T) f(x)\\r =\\ x - x ||} (66) 

and see the locality properties of the helds by plotting 
5^/(r) against r —\\ x — xq ||. For the anomaly part 
A(x), the variation of 7[ m)U! ] is considered. 

The results are shown in figures ^ El and 01 for P = 3.0. 
In figure ^ the variation of the topological held q(x) is 
shown. The locality of the topological held q(x) is clearly 
seen. We can read the locality range as g ~ 0.5. The 
maximum value of the held Aq(x) is also shown in the 
same hgure. We can conhrm that 

\Aq(x)\ ~O(10- 5 ) ~0(e- L / 2e ) (£ = 12). (67) 

In hgure|21 the variations of the current h^(x) are shown. 
It shows clearly that the current h fJl (x) has the same lo- 
cality property as the topological held q(x) has and thus 
is local. In hgure |3 the variations of the anomaly co- 
efficient J[ m ,w] an d the held Aq(x) are shown. We can 
conhrm that the gauge held dependence of j[ m ,w] is in- 
deed small and the same order of magnitude as the size 
of Aq(x) and its variation. 

The locality properties of the helds q(x), A(x), h^x) 
and Aq{x) are confirmed also for topologically non-trivial 
gauge helds. See hgure ^ 

We next examine the cancellation of the gauge 
anomaly. We consider the so-called 11112 model which 
consists of four Left-handed Weyl fermions with unit 
charge and one Right-handed Weyl fermion with charge 
two. The gauge anomaly cancellation condition in two- 
dimensions is satisfied as follows: 

4 

E e2 -( 2e ) 2 = - ( 68 ) 

In hgure we plot the anomaly parts of the Left- 
handed fermions 4 x A 1 (x) and of the Right-handed 
fermion A 2 (x), where A a (x) = A{x)\rj^ u<la , and the 
total anomaly part J2 a ^ a ( x ) — 4 x A (x) — A 2 (x) for 
(3 = 3.0 and m\2 = 0. The result is impressive. The 
size of the total anomaly part is reduced to the order 
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FIG. 1: The variation of the topological field 5 n q(r) (filled 
circle) is plotted against r =\\ x — xo ||. The maximum value 
of |Ag(x)| (cross) is also shown. The lattice size is L = 12. 
The gauge field is generated at /? = 3.0 with the vanishing 
magnetic flux ran = 0. 



FIG. 3: The variations 5 r ,7[ mjKI ] (filled diamond) and 5 v Aq(r) 
(filled triangle) are plotted against r =|| x — xo || . For a guide, 
the variation 5 v q(r) (open circle) and the maximum value of 
|A(jr(x)|(cross) are also shown. 
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FIG. 2: The variations S v hfj,(r) axe plotted against r —\\ 
x — xo || (filled square and triangle) . For a guide, the variation 
8 v q(r) (open circle) and the maximum value of |Ag(x)| (cross) 
are also shown. 



FIG. 4: The variations S v h IJ- (r) are plotted against r —\\ 
x — xo || (filled square and triangle). For a guide, the variation 
5 v q{r) (open circle) and the maximum value of |Ag(x)| (cross) 
are also shown. The lattice size is L = 10. The gauge field is 
generated at j3 = 3.0 with the magnetic flux mi2 = 1. 
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FIG. 5: The anomaly parts of the Left-handed fermions 4 x 
A 1 (x) (open circle) and of the Right-handed fermion A 2 (x) 
(cross) are plotted against the fused coordinate x + L x t. 
The total anomaly part 4 x A 1 (x) — A 2 (x) (filled rectangle) 
is also plotted. The lattice size is L = 10. The gauge field is 
generated at j3 = 3.0 with the vanishing magnetic flux m\2 = 
0. 



0(e L / 2e ) after the cancellation and we can confirm that 

AA{x) ~ 0(e~ L/2e ) ~ 0(lQ- 4 ) (L = 10). (69) 

This result is also confirmed in figure EJ In this figure, 
the total topological field ^2 a q a {x), the total anomaly 
part ^2 a A a (x) and the total finite volume correction 
^ A? a (a;) are plotted. We see clearly that the size 
of AA(x) is the same order of magnitude as the size 
of Aq(x). The similar cancellation is also observed for 
topologically non-trivial gauge fields with mi2 = 1. 

Now the current k^(x) can be computed from 
Eqs. (|55p. (|56|l and (|58|) . The current so obtained is 
thus local. This local current can reproduce the origi- 
nal topological charge q(x) within the deviation of order 
O(10 -15 ). The gauge invariance of the current is also 
maintained within the error. 



VI. DISCUSSION 

We have demonstrated that the cohomological anal- 
ysis of the chiral anomaly associated with the over- 
lap fermions can be performed numerically in two- 
dimensional lattice with a finite volume. The resulted 
current k^(x) is gauge invariant and local. 



FIG. 6: The total topological field J2 a <f{x) ( nll ed triangle), 
the total anomaly part ~}2 a A a (x) (filled rectangle) and the 
total finite volume correction ~}2 a Aq a (x) (open circle) are 
plotted against the fused coordinate x + L x t. 



Four-dimensional case is numerically demanding. The 
evaluation of the differential of the topological fields 
should be performed without diagonalizing H w . For 
other parts, our procedure would work also in this case. 

Our next step would be the numerical construc- 
tion of the Weyl fermions measures and observables 
in the anomaly- free U(l) chiral gauge theories in two- 
dimensions, keeping the gauge invariance. Work in this 
direction is in progress. 

APPENDIX A: 

In this appendix, we give the explicit formulae for the 
vector-potential representation of the link variables of the 
admissible U(l) gauge field discussed in the section ITT1 
We first define 

a^(x) = - ln{U(x, fi)V [m] (x, fi)^ 1 } (Al) 

and 

2nn^(x) = i^„(x) - 27T7 ^ U - {d^a v {x) - d v a^(x)}. 

(A2) 

We also introduce a summation convention as 

( eS/w (*i>i) 

E '/(*) = { ( Xi =0) . (A3) 

{ E^i x< (-i)/W 0*<-i) 
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Then, in four-dimensions, the vector-potential is de- 
fined by 



' X\ — 1 



X 2 — 1 



times -ff w is approximated by 



t 2 =o 



Ai(x) = a 4 (a;) + 27r< ^ 'nu{x)+ ^ '2nh 24 ,(x) | Xl=0 

Ui=0 

+ '™34(z)L=x 2 =0 f , 



/#2 



w fc=l 



hi + c 2 fe_i ' 



(Bl) 



t 3 =o 



' x\ — l 



A 3 (x) = h 3 (x) + 2irl ^ 'ni3(x)+ ^ '"23(a)|xi=o 

1*1=0 t 2 =0 

x 4 -l L/2-1 "I 

-^ 3 ,L/2-l Y ' Y "34(201x1=3:2=0 > , 
t i= t 3 =-L/2 J 

{Xl-1 
Y '"12( X ) 
ii=0 

x 3 -l L/2-1 

-4 2 ,l/2-i Y ' Y n 23 (x)\ Xl= o, Xi= o 

t 3 =0 t 2 =-L/2 

x 4 -l L/2-1 ] 

-£x 2 ,L/2-l Y ' Y "24(2;) |xi=0 > , 
t 4 =0 t 2 =-L/2 J 

Ai(a;) = ai(or) + 27r • 

x 2 -l L/2-1 

-^xi,L/2-l J! ' Y "12(2;)|x3=X4=0 
t 2 =0 ti=-L/2 
x 3 -l L/2-1 

-4i,L/2-l J! ' Y ^13(^)1x4=0 
t 3 =0 ti=-L/2 

x 4 -l L/2-1 

-4i,L/2-l J! ' X " 14 ( X ) 
t 4 =0 ti=-L/2 

In two-dimensions, the vector potential is defined by 

Xl— 1 

A 2 (x) = a 2 (x) + 2ir ^ 'n 12 (x), 

ti=0 

= ai(x) 

x 2 -l L/2-1 

-^xi,L/2-l X ' Y 2 ^i 2 {x). 
t 2 =0 ti=-L/2 

(A4) 



APPENDIX B: RATIONAL APPROXIMATION 
OF THE OVERLAP DIRAC OPERATOR 

In this appendix, we give the formula of the rational 
approximation of the overlap Dirac operator with the 
Zolotarev optimization. The inverse square root of 



h w is defined by H w /X m i D where A m j n is the square root 
of the minimum of the eigenvalues of H^, and the coeffi- 
cients Cfe and bk are given as follows: 



sm 



Cfe 



A = 



i-sn 2 (fr;0' 



d 



-TT — 
1 11 1 



fc=i 



C 2 k 



bk = d 



I\"=l( C 2t - C 2 fc-l) 
nr=l, 4 #*:( C 2*-l - c 2fc-l) 



(B2) 

(B3) 

(B4) 
(B5) 



K' is the complete elliptic integral of the first kind with 
modulus k' = y/l - (A min /A max ) 2 . 



APPENDIX C: SOLUTIONS OF THE MODIFIED 
POINCARE LEMMA AND THE COROLLARY OF 
DE RHAM THEOREM IN TWO-DIMENSIONS 



In this appendix, we give the explicit solutions of the 
Poincare lemma in two-dimensions (n=2). 
IV.a(O-form) and IV.b(O-form): 

Xl 

g 1 (x 1 ,x 2 ) = 5 X2iX o Y /(yi)' 

j/i=x°-L/2 

X 2 

g 2 {x u x 2 ) = Y {f( x uy2)-5 V2iX of(xi)}, 

y 2 =x° 2 -L/2 

(CI) 

where f(xi) = ^=-1/2 /0*i>i&)- 
IV.a(l-form): 

x 2 

9i 2 (xi,x 2 ) = - Y [h(xi,y 2 ) - ^ 2;X o/i(a;i)| , 

y 2 =x° 2 -L/2 
Af 1 (x 1 ,X 2 ) = ^x 2 ,x°/l(2;i), 

Af 2 (x!,X 2 ) = f2(xi,X 2 )\ X2 = x o_ L/2 -l, (C2) 

where /i(a;i) = E^Il/2 h{x\, y 2 )- 
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IV.b(l-form): 

9i2(xi,x 2 ) = - { h {xi , y 2 ■)- 6^^0/1(0:1)} 

y 2 =x°~L/2 

+ X] /2(yi,^2)U 2 = 2; 0-L/2-l- 
yi =lD°-L/2 

(C3) 
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